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Abstract. DENIS observations in the J (1.25 /im) and Ks 
(2.15 /im) bands together with isochrones calculated for 
the RGB and AGB phase are used to draw an extinction 
map of the inner Galactic Bulge. The uncertainty in this 
method is mainly limited by the optical depth of the Bulge 
itself. A comparison with fields of known extinction shows 
a very good agreement. We present an extinction map for 
the inner Galactic Bulge (~ 20 deg 2 ) 
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1. Introduction 

Studying the stellar populations in the Galactic Bulge re 
quires knowledge of the interstellar extinction. Previous 
studies (e. g. Catchp ole et al 
Unavane et al 



1990, Frogel et al. 1999 



1998) have shown that it can vary from 
Ay = 5 m up to Ay = 35 m towards the Galactic Centre 
region. As in most parts of the Galactic Bulge, interstel- 
lar absorption is not homogeneous but occurs in clumps, a 
detailed extinction map of the Bulge is therefore essential. 

Catchpole et al. ( 1990[ ) mapped the interstellar ex- 
tinction around the Galactic Centre (~ 2 deg 2 ) using the 
red giant branch of 47Tuc as a reference. Stanek et al. 
( |1996| ) mapped the interstellar extinction of Baade's win- 
dow using OGLE photometry of red clump stars. Frogel 
et al. ( |1999| ) determined the interstellar reddening for a 
few fields in the inner Galactic Bulge using the red gi- 
ant branch of Baade's window as a reference. Ak values 
ranging from 0.27 up to 2.15 mag were fo und b y the latter. 

The DENIS survey (Epchtein et al. |l997| , Persi et al. 
1995| ) with 2MASS (Skrutskie et al. |l997|) is the first at- 



tempt to carry out a complete survey of the southern sky. 

Send offprint requests to: schulthe@iap.fr 

* Based on observations collected at the European Southern 
Observatory, La Silla Chile 



The limiting magnitudes in the three near-IR photomet- 
ric bands (I = 0.8 /im, J = 1.25 /im and Ks = 2.15 /im) 
for point sources are 18 m , 16 m and 13 m respectively. The 
photometric accuracy (rms) is better than 0.1 mag and the 
astrometric accuracy better than 1 arcsec. These numbers 
are for uncrowded fields. 

DENIS Ks/(J — Ks) colour magnitude diagrams 
(CMDs) in regions of low extinction show a well defined 
RGB and AGB sequence in the Galactic Bulge (see Fig. 1). 
Due to their high luminosities (up to Mboi ~ —4.5 for 
non-saturated sources in DENIS) these stars are ideal 
tracers of the stellar populations in the Bulge and are 
found even in highly extincted regions. In this paper we 
present a method to derive interstellar extinction using 



isochrones from Bertelli et al. (1994) in combination with 
DENIS CMDs. We show that this method is appropriate 
for low as well as for highly extincted regions. Finally we 
present a map of the extinction in the inner Bulge between 
-8<^<8 and -1.5<b<1.5°. The finer details of the features 
seen in the map will be discussed in a subsequent paper. 

2. Observations 

The near infrared data were acquired in the framework 
of the DENIS survey, in a dedicated observation of a 
large Bulge field, simultaneously (Summer 1998) in the 
three usual DENIS bands, Gunn-I (0.8 /mi), J (1.25 /im) 
and Ks (2.15 /im). For the source extraction we used 
PSF fitting optimised for the crowded fields (Alard et al. 
in preparation). For the astrometry, the individual DE- 
NIS frames were cross-correlated with the PMM catalog 
(USNO-A2.0). The absolute astrometry is then fixed by 
the accuracy of this catalog (<~ 1"). The internal accuracy 
of the DENIS astrometry, derived from the identifications 
in the overlaps is of the order of 0.5" . 

For the determination of the zero point all standard 
stars observed in a given night have been used. The typical 
uncertainty of the zero points has been derived from the 
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overlapping regions and is about 0.05, 0.15 and 0.15 mag 
in the I, J and Kg bands respectively. 

3. Extinction determination using isochrones 

Colour-Magnitude Diagrams were constructed for sources 
in a small window (radius of 2 arcmin) in the field. 
The modal value of the distribution of the Ay re- 
quired to move the stars in the CMD to the zero ex- 
tinction isochrone was taken as the value of the ex- 
tinction for this window. The interstellar extinction law 
(A v : Aj : A Ks = 1 : 0.256 : 0.089) from Glass ( p.999| ) was 
used. The window was then displaced laterally in uniform 
steps to construct an image of the spatial distribution of 
the extinction over the whole field. We have found that 
towards the Galactic Bulge a sampling window of radius 
2 arcmin provides a sufficient number of stars to form 
a sequence enabling a reliable estimate of the extinction 
in such a window. We presently do not use the DENIS 
data with Kg < 7 due to the saturation of the detectors. 
We only use those sources which have been detected in 
J as well as in Kg with Kg brighter than 11 mag in or- 
der to be as complete as possible in J. This criterion is 
also quite important in order to rule out fake sources at 
the fainter end of the luminosity function and to guaran- 
tee an RGB/AGB identification. However, in the regions 
with very high extinction (Ay > 25) we find that a large 
proportion of the sources detected at Kg do not have coun- 
terparts at the shorter wavelengths. The location of these 
'missing' J sources is concentrated in regions with high 
extinction (Ay > 25 m ). Hence for regions with large ex- 
tinction, where the number of sources detected with the 
above criterion is smaller than in regions with lower ex- 
tinction, we only get a lower limit on the extinction. 

Results with such a sampling window and with dis- 
placement steps of 1 arcmin are discussed below. 

Isochrones from Bertelli et al. (1994) placed at 8kpc 



distance for a lOGyr stellar population with Z=0.02 has 
been used as a reference system. The isochrones were cal- 
culated for the ESO filter system by convolving the near- 
infrared bands with the spectral energy distributions from 
Kurucz ( 1992 ) for temperatures higher than 4000 K. At 
lower temperat ures the effective temperature scale from 
Ridgway et al. ( 1980 ) for the late M giants and the Lancon 
& Rocca-Volmerange (1992) scale for the early M giants 
has been used. The lack of very red standards limi ts the 
near-infrared colour transformation (Bressan & Nasi 1999 ) 
and causes the colours of the Z=0.02 isochrone to 'satu- 
rate' around (J — K)o = 1.35. Therefore a new, empirical 
Teff-(J-K) colour relation has been derived by making a 
fit through the T e ff/ (J — K )q data available for cool giants 
[see Schultheis et al. 1998| and Ng et al. (in preparation)]. 
Schultheis et al. (1998) demonstrated the good agreement 
of the isochrones with the new T e ff-colour relation using 
NIR photometry of a sample of Miras and Semiregular 
Variables in a field located at the outer Bulge. However, 
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Fig. 1. DENIS K S /(J - K) s diagram for Baade's win- 
dow (Sgr I). Superimposed are isochrones with Z=0.008, 
Z=0.02 and Z=0.05, respectively (see text). Saturated 
stars (Kg < 7) are not displayed. 



the upper part above the RGB tip (Kg ~ 8.0) remains 
nevertheless more uncertain. 

Based on observed near infrared spectra for a sam- 
ple of M giants and Mir a Variables, kindly provided by 
A. Lancon, we have found the difference between K and 
Kg to be small, in the order of 0.04-0.05 mag. 

3.1. Uncertainties of the isochrone method 

Figure shows a DENIS CMD in a part of the Baade's 
Window (SgrI, ~ 0.06 deg 2 ). Overlaid are the isochrones 
with different metallicities (for Ay = 0). The isochrone 
with Z=0.02 is found to follow well the observed CMD 
and is shown at Ay = and also at Ay = 1.35 which is 
the modal value of the extinction found for the field. 

It is a well known fact that in the Galactic Bulge 
one deals with a wide metallicity and age range. Using 
isochrones with a different Z would give us different extinc- 
tion values. In order to estimate the effect of metallicity 
on the derived extinction values, isochrones with Z=0.05 
and Z=0.008 have been considered. Between the isochrone 
with Z=0.008 and Z=0.05 there is a A(J - K) - 0.2 (see 
also e. g. Bressan et al. 1998) which corresponds to a 
AAy ~ l m . This is the typical uncertainty in the de- 
termination of the extinction assuming solar metallicity. 
In contrast to the effect of metallicity the isochrones are 
hardly affected by an age-range (see e. g. Schultheis et al. 
1998). 

Analysis of repeated observations (1996 & 1998) shows 
that the internal dispersion in the photometry, in the 
crowded regions, is less than 0.15mag for Kg < 11. 
Compared to the distance spread of the Galactic Bulge 
(~ 0.35 mag) the errors in the photometric accuracy as 
well as the errors coming from the isochrones are negligi- 
ble. 
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Fig. 2. Observed CMDs with isochrones at derived extinction values (left panel) and histograms of the extinction 
values for individual sources (right panel) for sample fields (radius=2') in Baade's window (filled circles, solid line), a 
field located at £=-0.37 and b — 0.5 (filled triangles, dotted line) and around the Galactic Centre (open circles, dashed 
line). 
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Fig. 3. Map of the extinction in the inner Galactic Bulge. The image grey scale represents an Ay range in 

magni tudes from (white) to 35 (black). A coloured high resolution image can be found at http: / /www 
denis.iap.fr/articles/extinction/ 



3.2. CMDs in sampling windows 

In Figure || we show the CMDs for sampling windows at 
three different locations, namely a field with low extinc- 
tion in Baade's window, one with intermediate extinction 
at £ = —0.37° and b = 0.5° and one at the Galactic 
Centre (£ = 0,6 = 0). The low extinction field consid- 
ered here is a small part of the Sgr I field whose CMD 
is presented in Figure 0. While the RGB/AGB branch 
in Baade's window shows a narrow RGB/AGB sequence, 
for the two other fields one sees clearly a wide spread in 
the (J — Kg) colour. The Baade's window is known to have 
low and well behaved extinction and this is seen as a sharp 
peak in the distribution of Ay (figure [2]) . For the field at 
£ — —0.37, b — 0.5 and around the Galactic Centre one 
does not find a single well defined peak but perhaps two 
different peaks which indicates that there maybe two or 
more distinct layers of extinction causing material along 
this line of sight. Alternatively the absorbing matter may 
show dumpiness on a scale smaller than the 4 arcminute 



diameter of the sampling window. As mentioned earlier, 
our J and K detections are complete in regions with low 
extinction (up to Ay < 25) and a large proportion of K 
sources do not have J counterparts (up to 60-90%) in the 
regions with high extinction (GC). It is unfortunately not 
possible to use the results from the regions with low ex- 
tinction to assign a completeness limit for the obscured 
ones because the controlling factors are confusion in the 
first case and the detector sensitivity in the second. A map 
of the 'missing' J sources shows that they are a significant 
source of uncertainty in regions with Ay > 25. 

4. Results 

The whole map for the inner Galactic Bulge is shown in 
Fig. | The map has a resolution of 4 arcminutes (the 
diameter of the sampling window used). Note the clumpy 
and filamentary behaviour of the distribution of extinction 
especially close to the Galactic plane. There are also small 
pockets with uniform extinction of about Ay = 6 magni- 
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Table 1. Mean value of A v (in 0.1° x 0.1° box) along 
different lines of sight towards the inner Galaxy 
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tudes (for example the ISOGAL field I = 0, b = 1 studied 
by Omont et al. 1999). Table | gives the values of the mean 
extinction in various locations in the inner Bulge. The dis- 
tribution of extinction along the line of sight through the 
Bulge will be discussed in detail in a subsequent paper 
(Ganesh et al in preparation). 

In Fig. ^ we present a contour map of the extinction 
around the Galactic Centre. The distribution of extinc- 
tion in this region has been studied earlier by Catchpole 
et al ( 1990D using H&K data. We find a good agreement in 
the observed structures keeping in mind the higher resolu- 
tion results presented here. The difference in Ay compared 
to Catchpole et al. is typically smaller than 3. The aver- 
age interstellar extinction along the Galactic plane (see 
Fig. Q is w 24 m in Ay although this might be a lower 
limit (as seen by the presence of patches with Ay > 30) 
due to the effect of the 'missing' J sources discussed ear- 
lier. Along the minor axis the Ay decreases to 15 m towards 
the edges, starting from nearly 25 m near the Centre. The 
derived extinction values from Fig. |J are also qualitatively 
in agreement with the values found by Wood et al. (1998) 
near the Galactic Centre. 



5. Conclusion 

DENIS observations in the J and Ks band together with 
isochrones by Bertelli et al. (1994) are an excellent tool 



to map the interstellar extinction. A comparison with the 
field studied by Catchpole shows a very good agreement 
although we see more details due to the better resolution. 
Several fields with relatively low (Ay ~ 6) and homo- 
geneous extinction can be identified on the basis of the 
extinction map. This identification should help further de- 
tailed investigation of the stellar population in these win- 
dows. A study of the three-dimensional distribution of the 
material responsible for the interstellar extinction should 
be facilitated by the availability of the extinction map. 
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Fig. 4. Contour plot of the extinction around the Galactic 
Centre with contours plotted in the range from to 35 at 
7 mag intervals in Ay. The major contours are labeled in 
large italic fonts. Local maxima (with Ay > 30) are also 
marked with the local Ay values in small fonts. 
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